http://www.arocmag.com/article/02-2019-10-008.html 


基于 非 线性 渐变 式 颜 色 映射 的 LIC 改进 算法 “ 


(上 海 海事 大 学 信息 工程 学 院 ， 上 海 201306) 


摘 要 : ARTAR (LIC) 算法 展现 的 纹理 反映 了 整个 失 量 场 的 方向 结构 ， 但 却 不 能 展现 失 量 场 的 强度 大 小 。 针 对 
此 问题 提出 基于 非 线 性 渐变 式 颜 色 上 映射 的 LIC 改进 算法 ， 将 矢量 场 强度 与 白 骂 声 结合 作为 LIC 的 输入 纹理 ， 运 用 
FastLIC 思想 并 划分 纹理 区 域 同步 执行 LIC 运算 来 提高 算法 效率 ; 再 将 失 量 场 强度 作 非 线性 变换 ， 根 据 渐 变 式 颜 色 映 
射 方案 使 用 OpenCV 处 理 引擎 并 行 实现 矢量 场 强度 的 颜色 映射 ; 最 后 由 LIC 得 到 的 灰 度 纹理 和 颜色 映射 结果 确定 合 
成 系数 并 构造 累计 函数 增强 二 者 结果 ， 再 进行 线性 合成 运算 得 到 最 终 的 可 视 化 效果 。 在 对 全 球 海洋 流 场 和 风 场 两 种 典 
型 的 矢量 场 进行 可 视 化 以 及 与 其 他 算法 对 比 实验 表明 ， 改 进 算 法 得 到 的 可 视 化 效果 纹理 清晰 ， 较 好 地 展示 出 矢量 场 的 
方向 和 强度 ， 能 够 更 准确 地 反映 和 失 量 场 的 全 方位 信息 和 局 部 变化 情况 。 
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Improved LIC algorithm based on nonlinear gradual-changing color mapping 


Gao Maoting, Dong Hongzhou, Zhou Fan 
(College of Information Engineering, Shanghai Maritime University, Shanghai 201306, China) 


Abstract: The texture revealed by line integral convolution (LIC) reflects the directional structure of the whole vector field, 


but cannot show the field intensity. In order to solve this problem, this paper proposed an improved LIC algorithm based on 

© nonlinear gradual-changing color mapping. Firstly, the intensity of vector field combined with white noise to form the input 

C texture of improved LIC, which used the idea of FastLIC and divides the texture into several regions for synchronous 
execution of LIC, to improve the efficiency of LIC. Secondly, the improved LIC does a non-linear transformation on the 
intensity of vector field, and used the processing engine, OpenCV, to implement the color mapping of the vector field intensity 
according to the gradual-changing color-mapping scheme. Finally, the results of gray texture obtained by LIC and color 
mapping determine the synthetic coefficient, constructing a cumulative function enhances the two results, and adopting the 
way of linear synthesis obtains the final visualization. Simulation results show that when applied to the global ocean flow field 
and wind field, the proposed algorithm generates a clear visualization displays the direction and intensity of vector field better, 
and reflects the global information and local changes of the vector field better, when compared with other algorithms. 
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0 ala 场 的 方向 变化 和 形状 特征 ， 在 矢量 场 可 视 化 中 具有 独特 的 优势 
e w, 

采用 可 视 化 的 方式 对 矢量 场 数据 展示 往往 是 掌握 和 分 析 其 对 矢量 场 数 据 可 视 化 需要 同时 考虑 方向 和 强度 大 小 ， 但 

运动 规律 的 重要 手段 ， 也 为 进一步 发 现 矢量 场 的 特征 结构 提供 LIC 算法 生成 的 纹理 只 能 展示 矢量 场 的 方向 ， 并 不 能 反映 出 矢 


了 基础 支撑 。 但 矢量 数据 因 具 有 大 小 和 方向 两 个 属性 ， 目 前 还 量 场 的 强度 信息 ， 而 且 算法 比较 耗 时 。 针 对 LIC 算法 存在 的 问 
没有 直接 适用 的 可 视 化 模型 叫 。 实 际 应 用 中 的 矢量 场 ( 如 海洋 流 。 B, 许多 改进 算法 先后 被 提出 ， ue 了 基于 HLS 模 


场 、 风 场 等 ) 往 往 具 有 复杂 性 和 特殊 性 ， 传 统 的 点 图 标 法 和 矢量 “型 的 FastLIC 算法 ,将 矢量 大 小 通过 HLS 颜色 模型 映射 ， 再 利 
线 法 中 已 经 无 法 满足 其 可 视 化 要 求 。 而 线 积分 卷 积 算法 Cine 。 用 OpenGL 技术 将 映射 结果 三 加 到 FastLIC 算法 生成 的 灰 度 纹 
integral convolution, LIC) 通过 生成 连续 平滑 的 纹理 展现 矢量 HE, X LIC 算法 展现 矢量 场 强度 提供 了 新 思路 ， 陆 剑 锋 等 人 
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k 


第 一 个 模块 是 LIC 运算 的 改进 ， 将 矢量 场 强度 与 白 噪声 相 由 于 和 矩阵 了 中 的 rand(x, y) e[0,1] ,为 了 和 白 噪 声 与 矢量 场 强 


结合 形成 新 的 纹理 作为 LIC 运算 的 输入 , 在 LIC 得 到 的 灰 度 纹 ” 度 更 好 地 融合 ， 此 处 应 该 把 场 强 和 矩阵 NN 中 的 值 归 一 化 到 [0 
里 中 融入 场 强 信息 进 行 增强 显示 。 为 了 提高 LIC 的 运算 速度 ， 内 ， 即 令 式 (2) 中 a=0,b=1。 由 式 (1) 得 LIC 算法 输出 的 纹理 
降低 算法 对 硬件 的 依赖 程度 ， 该 模块 在 FastLIC 算法 思想 的 基 强度 O) 与 白 噪 声 7 成 正 相 关 ， 在 融入 场 强大 小 后 ，7(x) 和 
出 上 ， 将 合成 纹理 分 成 若干 区 域 同步 执行 LIC 运算 。 新 生成 的 纹理 依然 要 保证 这 种 相关 性 。 所 以 本 文采 用 将 白 噪 声 

第 二 个 模块 是 非 线性 渐变 式 颜 色 映 射 ， 根 据 实际 应 用 场景 EREET 与 归 一 化 的 矢量 场 强度 矩阵 IN. 进行 加 权 求 和 的 方式 , E. 


选 定 颜色 表 ， 制 定 非 线性 渐变 式 颜色 映射 方案 并 采用 OpenCV — 阵 了 中 的 各 项 取 值 为 [0] ，M 归 一 化 后 各 项 数值 的 区 间 为 
实现 ， 克 服 了 线性 映射 在 矢量 场 强度 变化 不 明显 区 域 映 射 效果 [0 。 根 据 和 矩阵 7 ，XWN 定义 一 个 生成 新 的 输入 纹理 的 函数 


不 佳 的 缺点 。 该 模块 与 第 一 模块 将 实现 并 行 计算 来 提高 算法 的 。””g(7T,N), 令 权重 系数 为 &,B，, 白 噪 声 融入 场 强大 小 后 形成 的 新 
运行 速度 。 的 纹理 为 2 ， 则 

第 三 个 模块 是 灰 度 纹理 与 非 线 性 渐变 式 颜 色 映 射 结果 合 gopr We hy 加 
成 ， 首 先 根据 各 个 点 的 灰 度 纹理 值 和 映射 结果 确定 线性 合成 系 arp 
数 ， 再 通过 构造 附加 的 累计 函数 ， 将 灰 度 纹理 值 和 映射 结果 根 a 越 大 ， 表 示 白 噪声 所 占 比重 越 高 ， 由 此 生成 的 输出 纹 忆 
据 函 数 作 相应 变换 ， 最 后 合成 运算 得 到 最 终 的 彩色 纹理 结果 ， 将 越 能 体现 矢量 场 的 方向 信息 ;有 越 大 ， 表 示 矢 量 场 强 度 所 占 
克服 了 一 般 线 性 合成 方法 所 得 效果 亮度 暗 、 纹 理 不 清晰 以 及 合 比重 越 高 ， 最 终生 成 的 可 视 化 效果 越 能 体现 出 矢量 场 的 强度 信 
成 系数 难以 确定 等 缺点 。 息 。 
2.1 LIC 运算 的 改进 获得 新 的 输入 纹理 后 ，LIC 算法 需要 从 当前 点 开始 分 别 向 


经 典 LIC 算法 获取 原始 矢量 场 后 ， 生 成 一 个 与 原始 矢量 场 正 、 反 两 个 方向 延伸 到 两 段 长 为 工 的 流 线 ， 此 过 程 称 为 流 线 生 
规模 一 致 的 白 噪声 作为 输入 纹理 。 假 设 二 维 矢量 场 大 小 为 。” 成 (1。 本 文选 用 数值 积分 的 方式 生成 流 线 。 在 对 海洋 流 场 和 风 
和 x7 ， 由 原始 矢量 场 得 到 的 白 噪声 纹理 为 7， 为 了 计算 方便 ， ” 场 可 视 化 的 应 用 中 ， 由 于 实测 数据 是 基于 采样 点 的 ， 这 里 需要 
] rand(x. y) 表示 了 中 各 个 项 的 取 值 ，rand(x,y) 73r T [0.1] 2 用 离散 的 思想 描述 LIC 算法 中 流 线 的 轨迹 s(?) ， 用 V(s(D),) 表 


Sr 


间 的 随机 数 ， 则 可 令 T 为 示 流 线 速度 ， 则 在 一 个 时 间 段 [sh] 内 , 流 线 的 轨迹 可 表示 为 
rand(ll) rand(l,2) --- rand(ly) 0, = s) + [EVON i$ 
rand(2,)) rand(2,2) -.. rand(2, y) i 
"e s l s. i 其 中 : VGG.D 又 可 以 分 解 为 
V(Gs(t).t) 2 [v, (x, y],7),v, (o y],7)] (5) 


rand(x,l) rand(x,2) … rand(x, y) 


对 式 (4) 进行 求解 即 可 得 到 该 点 生成 的 流 线 。 为 了 对 矢量 场 
文献 [15] 提 出 将 白 噪 声 纹理 着 色 来 区 分 不 同 的 运动 方向 ， 中 复杂 微小 结构 进行 很 好 地 展示 ， 本 文选 用 精度 最 高 的 变 长 
本 文 方法 在 白 噪声 纹理 中 融入 原始 矢量 场 的 强度 信息 进行 LIC ffr Runge-Kutta 方法 08 求 出 生成 流 线 ， 而 且 随 着 矢量 场 数据 规 
运算 。 由 式 (1) 可 得 在 白 噪声 纹理 中 融入 场 强 大 小 后 形成 新 的 输 — 模 的 增加 ， 基 于 RK4 的 求解 方法 加 速效 果 明显 。 
入 纹理 完全 可 行 。 通 过 原始 矢量 场 获得 场 强 矩阵 后 ， 为 了 能 然而 流 线 生 成 模块 因 重 复 计 算 大 量 流 线 导 致 LIC 算法 十 分 
TE ELSE PRAE T 中 融入 场 强大 小 ， 需 要 对 原始 矢量 场 强度 矩阵 耗 时 。 为 了 提高 算法 的 运行 效率 , 本 文 对 LIC 运算 进行 了 改进 。 


进行 归 一 化 处 理 09。 定 义 原始 矢量 场 强 度 和 矩阵 为 M ， 归 一 化 前 降低 LIC 算法 时 间 复 杂 度 的 方式 主要 分 为 四 类 : 第 一 类 是 
处 理 后 的 算 阵 为 V ， 现 给 出 将 矩阵 M 中 所 有 项 归 一 化 [4,D] 区 ”将 LIC 算法 改 为 FastLIC 算法 ， 不 用 重复 计算 流 线 ， 而 是 充分 
间 的 一 般 方法 : 利用 了 点 之 间 的 相关 性 ， 即 后 一 点 的 纹理 值 可 以 由 前 一 个 点 得 
MET m 出 ; 第 二 类 是 并 行 LIC 算法 ， 在 大 规模 并 行 分 布 式 存储 计算 机 
au C Main 上 实现 ， 速 度 快 但 对 硬件 要 求 特别 高 ， 第 三 类 是 分 级 线 积分 策 
Kp. m, 为 矩阵 M 中 的 任意 项 ，7 为 ws 归 一 化 后 对 应 的 项 ， 。 略 吧 ， 利 用 多 核 架 构 对 算法 中 的 积分 过 程 进行 加 速 ， 缩 短 LIC 
ms Elab]; m, 为 矩阵 M 中 项 的 最 大 值 ，mam 为 矩阵 M 中 项 ”运算 时 间 ， 第 四 类 则 是 通过 GPU 实现 LIC 算法 的 并 行 计算 ， 
的 最 小 值 。 最 后 所 得 矩阵 N 为 利用 GPU 大 量 的 执行 单元 提高 算法 效率 。 第 二 类 和 第 四 类 对 硬 
"PR 件 要 求 过 高 ， 对 于 第 三 类 改进 算法 ， 如 果 矢 量 场 规模 增加 ， 则 
运算 时 间 依旧 呈 线 性 增长 。 本 文选 择 FastLIC 算法 思想 来 提高 
TT LIC 算法 的 效率 ， 并 在 此 基础 上 将 输入 纹理 划分 为 个 区 域 ， 
TOROS 分 模块 同步 执行 LIC 运算 。FastLIC 算法 思想 提出 ， 根 据 相 邻 
na no o My 位 置 点 之 间 的 相关 性 ， 可 以 通过 已 知 点 的 灰 度 值 求 其 相 邻 点 的 


灰 度 值 ， 避 免 不 必 要 的 流 线 生 成 计算 。 假 设 羡 和 浆 是 矢量 场 上 
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相 邻 的 两 点 ， 计 算 它 们 的 卷 积 时 ， 可 以 发 现 这 两 点 4 


有 很 大 部 分 相互 覆盖 的 区 域 ， 如 图 3 所 示 。 


过 式 (6) 求 得 


J3 ”点 和 x 点 生成 的 流 线 
若 已 知 处 的 灰 度 值 为 1(%) ， 则 交 处 的 灰 度 值 O) 可 通 


I(x) 2 I0) - IA * I(A;) 


式 (6) 只 在 盒 式 卷 积 核 的 情况 下 成 立 ， 


别 是 A 和 A， 


区 域 的 卷 积 结果 。1 


性 ， 通 过 式 (1) 可 得 到 采样 点 离散 情况 下 为 点 处 纹理 灰 


Léa) 的 计算 方法 如 式 (7) 所 示 。 


I(x)- X Tk) 


其 中 : x 为 流 线 上 离散 点 ; T) 73 x; AAAA AE, 


78 x X xs BOO FL. f ka) =L. 


点 离散 的 特 ' 


根据 采样 
则 正 反 流 线 上 共 


成 的 流 线 


(6) 
H (A) 7 ICA;) 分 
于 数据 采样 点 具有 离散 的 特 


度 值 


() 


k(x;) 


性 , 若 每 条 长 为 荆 的 流 线 上 有 并 个 点 ， 
J 2n+1 个 点 。 对 于 离散 点 的 卷 积 ， 


需要 满足 


所 有 参与 卷 积 的 点 的 卷 积 系数 之 和 为 1， 且 每 个 点 的 系数 都 相 


等 


所 以 得 到 系数 为 1/(22+JD ， 即 式 (7) 


FP 的 K(x)=1/ (2n+D),， 


满足 盒 式 卷 积 核 ， 不 仅 保 证 了 显示 质量 同时 得 到 最 快 的 运算 
E, FE x 点 处 的 灰 度 值 为 
Iwz ÈT (8) 
YES (QD FRU 138 on 点 处 的 灰 度 值 计算 为 
u OT (x) + BN(x) 
1) 7 7- no "T ) (9) 
所 以 ， 运 用 FastLIC 思想 ， 基 于 盒 式 卷 积 核 ， 对 于 任意 从 
点 和 开始 的 LIC 运算 ， 不 需要 重复 计算 大 量 流 线 ， 只 要 求 出 
Ix) 的 值 ，7GCo 和 7(%, 1) 的 值 便 可 通过 式 (10) (1D 求 得 ， 
且 求 解 采 用 骨 入 式 变 步 长 龙 格 塔 库 数值 积分 法 ， 提 高 了 流 线 计 
算 的 准确 性 。 


TX) 2 I,) +K Om) -QGo 1» 
ana) =L En) +K Ona) -QG4) 


式 (10) (11) 中: 


k=1/ Qn1) ; 


2.2 非 线性 渐变 式 颜 色 映射 
为 了 更 好 地 展现 


量 场 的 强度 信息 ， 
文通 过 自选 定 颜色 表 制 定 非 线 性 渐变 式 颜 色 映射 方案 对 和 量 场 
强度 进行 颜色 泻 染 并 采用 OpenCV EH. OpenCV 是 一 个 路 平 


HW 
ES 
a 
Md 
3E 
ra 
nm 
HH 


(10) 
(41) 


QG 为 式 (3) 计算 得 到 的 新 纹 


， 本 
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台 的 计算 机 视觉 库 , 其 强大 的 视觉 处 理 算法 能 够 准确 地 从 选 定 
的 颜色 表 中 获取 颜色 样本 ， 不 需要 添加 新 的 外 部 支持 也 可 以 完 
整地 编译 链接 生成 执行 程序 ， 在 颜色 映射 方面 处 理 十 分 高 效 。 


本 文 将 借助 OpenCV 轻 量 级 处 理 引 擎 快速 实现 非 线性 渐变 式 颜 


色 映 射 方案 为 了 进一步 提高 算法 的 效率 , 颜色 映射 过 程 与 LIC 
运算 采用 并 行 的 处 理 方式 。 


2.2.1 非 线 性 场 强 值 转换 

在 实际 应 用 中 ， 矢 量 场 强度 往往 分 布 不 均匀 ， 如 果 场 强身 
中 在 一 段 较 小 的 范围 内 ， 线 性 映射 会 将 场 强 映射 结果 集中 在 某 
一 种 颜色 区 域内 ， 这 样 就 很 难 区 分 出 局 部 场 强 的 大 小 。 针 对 线 
性 映射 存在 的 缺点 ， 本 文通 过 非 线 性 方式 转换 场 强 值 。 从 原始 
矢量 场 中 得 到 场 强 大 小 mag ， 将 场 强 大 小 根据 式 (12) 进行 转 
换 : 


A 


ns 1 "T 


newmag — 
c-1 

Hop. newmag 为 非 线 性 变换 后 的 结果 ; c 为 非 线性 映射 因子 ， 

实验 一 般 取 c=0.1; magw 为 矢量 场 强 度 中 的 最 大 值 。 根 据 

c=0.1 得 到 式 (12) 图 像 如 图 4 所 示 ， 由 图 像 可 得 


newmag € [0.1] 。 


Z] 4 c=0.1 If, newmag 与 mag 的 映射 关系 
2.2.2 渐变 式 颜 色 映射 

由 于 本 文 的 实验 对 象 是 海洋 流 场 和 风 场 ， 为 了 得 到 较 好 的 
颜色 映射 效果 , 本 文 的 颜色 模型 将 参照 PanoplyWin29 软 件 中 提 
供 的 颜色 表 。 in cdi 软件 是 由 NASA.GISS 发 布 的 海洋 数 
据 可 视 化 软件 ， 拥 有 专业 的 海洋 要 素颜 色 映 射 体系 。 从 
PanoplyWin 软件 中 选取 对 应 海洋 流 场 和 风 场 映射 的 颜色 表 如 
图 5 所 示 。 


E —m 


5  PanoplyWin 软件 中 对 应 流 场 和 风 场 映射 的 颜色 表 
在 原始 场 强 经 过 非 线性 变换 得 到 新 的 场 强 值 后 ， 为 了 实现 
颜色 映射 的 平滑 过 渡 ， 本 文 提出 如 下 的 渐变 式 颜色 映射 方案 : 


N 


OpenCV WREEK wxh, 若 要 将 非 线 性 变 
换 后 的 newmag 根据 颜色 表 进 行 映 射 ， 还 需 将 newmag 按照 式 
(13) 作 放大 处 理 : 


p -|[newmag x w | 

p jj newmag 变换 后 的 值 ，p e[0,w] 。 将 颜色 表 的 长 、 宽 

分 别 等 分 成 w、h 份 ， 取 其 长 度 上 任意 相 邻 的 两 点 ， 记 为 
Pm Pms B P E ps Pma Z, P 到 p, RESTI, P 


(13) 
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到 Pmt 的 距离 为 L , 如 图 6 所 示 。 


PEPEN MCI NN 
Ó—— —6$— — oO 
Pn p Pari 


图 6 在 颜色 表 的 长 度 上 任 取 两 个 相 邻 的 点 
若 要 达到 渐变 效果 ， 则 P 处 的 颜色 取 值 可 由 式 (14) (15) 得 


l l 
l =—— *ol +—— *col 
color| p] A color| p,,] LL color| Paa] (14) 


hl 
color[p,,]= Rotorp, ,可 (15) 


XX (14) (15) 中 : color p, i] 7326. E Cp, D) 处 的 颜色 值 。 将 4 
和 4 的 取 值 代入 式 (14) 得 到 原始 矢量 场 中 强度 mag 对 应 的 颜 
色 映 射 值 为 


color[ p] = (Pn — newmag * w)* color| p,,] 


(16) 


+ (newmag * w— p, colori p,a] 


2.3” 灰 度 纹理 与 非 线性 渐变 式 颜 色 映 射 结果 合成 
在 两 个 模块 分 别 执行 完 后 ， 需 要 进行 合成 才能 得 到 最 终 的 
可 视 化 效果 。 令 LIC 得 到 的 灰 度 纹理 为 Rijjc ， 颜 色 映 射 结果 为 
Rio, ， 常 见 的 合成 方式 如 式 (17) (18) 所 示 。 
Ros =tRic + (0 - DR, (17) 
=1x Rye X Ro (18) 


R, 


last 


式 (17) 08) F: Ru 为 合成 的 最 终结 果 ; 1 为 合成 系数 ， 且 
fs(0.D 。 

通过 式 (17) 的 合成 ， 虽 然 可 以 通过 参数 1 来 控制 灰 度 纹理 
和 颜色 映射 结果 的 显示 比重 ， 但 得 到 的 效果 整体 亮度 偏 暗 ， 对 
比 度 、 视 觉 分 辨 度 降 低 , 局 部 纹理 和 颜色 不 清晰 等 。 通过 式 (18) 


aXivA (ERAH 


Chi 
BOXER, S: ATIRAN as E RU ; 


假设 Ric 中 有 两 个 像素 点 Bs,B， 且 1(R,)<1(R,)， 
hp(L) ，i=123… 汉 在 定义 域内 为 增 函 数 ， 将 及 ,及 ,的 
入 上 述 式 (19) 中 , 得 7\(R,)<7'(R,)， 所 以 通过 附加 累计 函 
强 显示 的 过 程 并 不 改变 原来 灰 度 纹理 的 变化 趋势 。Rw: 中 
方式 同上 述 Ruc 的 方法 , 结果 记 为 上 (xm) , 则 最 终 的 合成 方 


I(x) Ij I.(x) 
T(x0) * LG) IQ) t LO) 


其 中 : LO) 为 非 线性 渐变 式 颜 色 映射 后 的 结果 值 。 
3 ”实验 结果 及 其 分 析 
3.1. 数据 集 和 实验 环境 


T's (X0) = 


因为 
ia 
数 增 
处 理 
式 为 


L'o) (20) 


本 文选 取 全 球 海洋 流 场 和 风 场 这 两 个 典型 的 矢量 场 作 
念 对象， 采用 由 美国 国家 海洋 和 大 气管 理 局 (NO A A) 提 人 4 
球 海洋 流 场 和 风 场 实测 数据 。 数 据 格式 为 NetCDF， 时 间 


的 全 


分 辨 


率 为 6 h， 空 间 分 辨 率 为 0.25°x0.25* F A 26 v 


EE] 
H AN 


180° W~180° E, 90° S~90° N. 


本 文 的 实验 环境 为 Windows 10 操作 系统 ，Intel core i5 处 


里 器 和 8 GB 内 存 , 开发 工具 为 VS2013, 图 形 
代码 使 用 C++ 语言 实现 。 

3.2 可 视 化 效果 对 比 实验 

本 文选 用 适合 海洋 流 场 、 风 场 的 专业 颜色 映射 表 得 到 
视 化 效果 与 基于 HLS, HSV 等 模型 映射 的 LIC 算法 得 到 的 
化 效果 存在 较 大 的 颜色 差 ， 区 分 度 很 大 。 所 以 本 实验 是 在 
颜色 映射 表 下 进行 ， 将 本 文 算 法 三 个 模块 得 到 的 可 视 化 结 
别 同 对 应 的 经 典 LIC 算法 、 线 性 颜色 映射 方法 以 及 基于 线 | 
色 映 射 的 LIC 算法 (LCLIC) 所 得 结果 进行 比较 ， 以 此 来 验 i 


库 选 用 OpenCV, 


的 可 
可 视 
同一 
果 分 
生 颜 
证 本 


文 算法 在 矢量 场 可 视 化 应 用 中 的 优 劣 性 。 选 用 全 球 海洋 流 ] 


a 


的 合成 ， 不 仅 同 时 弱化 了 纹理 和 颜色 映射 的 效果 ， 而 且 运算 速 
度 大 大 降低 。 由 此 可 知 ， 上 述 两 种 方式 的 合成 效果 并 不 理想 。 

本 文通 过 构造 附加 累计 分 布 函数 来 增强 灰 度 纹理 与 颜色 映 
射 结果 的 显示 。 过 程 如 下 : 

a) CI Ruc 的 大 小 为 X*Y， 得 到 Ric 中 的 最 大 灰 度 值 为 
lycsas ， 最 小 灰 度 值 为 Tourcnin o 

b) 将 Re 中 的 灰 度 划分 为 n 个 等 份 ， 并 标号 为 
LL Los L, ^ SEA DX HS REA Guess -Tucmn)/n。 

c) 通 过 统计 获得 Ric FUKEB YE 4E BL, Loss L, PAR 
个 数 ， 记 为 有 L),i= 上 ,2,3,…,n ， 并 分 别 求 出 和 L) 占 总 数 的 百 
DHE, WA AL) =AL)/(X *Y), i=1,2,3, n. 

d) 计算 各 个 灰 度 级 的 累计 分 布 aL), B 


hp(L) - 3 hs(L), i=1,2,3, 5n 
j=0 


e) 重新 分 配 灰 度 值 x) ， 因 为 灰 度 值 的 取 值 为 [0,255] , 
所 以 实验 时 nn 取 255， 即 


iL, 
2555 URSI MN JN 

7 (xzo) = yo x*y (19) 
0 L-0 


J 


可 视 化 过 程 进 行 实验 说 明 ， 最 后 再 给 出 本 文 算 法 应 用 到 全 


球 海 


洋 风 场 中 得 到 的 可 视 化 效果 。 实 验 采 用 控制 可 变 因素 的 方 
行 ， 经 过 前 期 多 次 实验 比较 ， 当 划分 模块 的 个 数 n=5 时 ， 
较 好 地 平衡 处 理 速 率 和 可 视 化 效果 质量 ， 积 分 半 长 L= 
a-B-M2 时 ， 可 视 化 效果 最 佳 。 所 以 ， 对 比 实 验 将 以 此 为 
进行 。 
在 对 全 球 海洋 流 场 数据 可 视 化 的 应 用 中 ， 本 文 改 进 的 
运算 模块 输出 的 结果 如 图 7(a) 所 示 , 经 典 LIC 算法 输出 的 结 


式 进 
能 够 
10, 

基准 


LIC 


如 图 7(b) 所 示 ， 黑 色 阴 影 部 分 为 陆地 面积 。 本 文 改进 的 LIC 38 


算 模 块 得 到 的 可 视 化 效果 ， 已 经 通过 明暗 渐变 的 纹理 初步 显 
下 也 十 分 清晰 。 


显示 


Ca) 本 文 改进 的 LIC 运算 模块 得 到 的 灰 度 纹理 效果 


录用 稿 高 茂 庭 ， 等 : 基于 非 线性 渐变 式 颜 色 映 射 的 LIC 改进 算法 


(b) 经 典 LIC 算法 得 到 的 灰 度 纹理 效果 

图 7 ”本文 改进 的 LIC 运算 模块 与 经 典 LIC 算法 的 灰 度 纹理 对 比 

图 8(a) 为 本 文 非 线 性 渐变 式 颜色 映射 的 效果 ; b) 为 采用 
线性 方式 颜色 映射 的 效果 。 本 文 提 出 的 非 线性 渐变 式 颜色 映射 
方法 能 够 适应 局 部 流 场 强度 变化 不 明显 的 区 域 。 当 局 部 区 域 流 
场 强度 变化 缓慢 时 ， 普 通 的 线性 映射 方法 会 在 该 区 域 形 成 较 大 
一 片 颜 色相 同 的 区 域 ， 如 图 8 (b) 中 的 方 框 标记 区 域 ， 而 基于 本 
文 的 方法 则 可 以 反映 出 流 场 强度 变化 较 小 的 分 布 情况 ， 如 图 
8 (a) 方 框 标 记 区 域 。 


(b) £207 时 LCLIC 得 到 的 全 球 流 场 可 视 化 效果 
图 9 线性 合成 不 同 合成 系数 对 应 的 可 视 化 效果 


Ca) 非 线性 渐变 式 颜色 映射 


图 10 本文 算法 得 到 的 全 球 流 场 可 视 化 效果 
海洋 风 场 可 视 化 的 处 理 流程 同 海洋 流 场 ,图 11 为 本 文 算法 
得 到 的 全 球 海洋 风 场 可 视 化 效果 。 通 过 海洋 流 场 和 风 场 这 两 个 
复杂 矢量 场 的 可 视 化 效果 来 看 ， 本 文 提 出 的 基于 非 线性 渐变 式 
颜色 映射 的 LIC 改进 算法 ， 运 用 在 矢量 场 可 视 化 中 ， 不 仅 清晰 
地 展现 了 矢量 场 的 方向 变化 ， 还 能 表现 出 矢量 场 的 强度 信息 。 


可 视 化 效果 具有 较 高 的 亮度 和 对 比 度 ， 克 服 了 其 他 LIC 改进 算 
(b) 普通 的 线性 映射 法 得 到 的 彩色 纹理 饱和 度 下 降 、 色 彩 偏 暗 、 纹 理 区 分 度 较 低 等 
图 8 非 线性 渐变 式 颜色 映射 与 普通 线性 映射 效果 对 比 缺点 ， 在 矢量 场 局 部 区 域 的 可 视 化 中 ， 也 能 够 清晰 地 展现 出 方 


LCLIC 算法 是 采用 式 (17) 的 合成 方式 , 在 这 种 合成 方式 中 ， 向 和 强度 的 变化 趋势 。 
合成 系数 并 没有 明确 的 规定 方法 。 图 9(a) 为 式 (17) 中 上 = 0.3 时 
LCLIC 算法 得 到 的 可 视 化 结果 , 纹理 不 清晰 , 但 色彩 比较 丰富 ; 
图 9(b) 为 式 (17) P £ 20:7 时 LCLIC 算法 得 到 的 可 视 化 结果 , ex 
理 效果 比 颜 色 效 果 明 显 。 除了 合成 系数 影响 外 ,基于 LCLIC 算 
法 得 到 的 可 视 化 效果 仍然 会 整体 变 暗 ， 纹 理 和 颜色 在 合成 后 效 
果 都 会 被 减弱 。 本 文 算法 明确 了 合成 系数 ， 并 构造 附加 累计 函 
数 对 两 者 结果 进行 变换 ， 克 服 了 合成 系数 难以 确定 以 及 合成 效 
果 不 理想 等 缺点 。 图 10 为 本 文 算法 得 到 的 全 球 流 场 可 视 化 效 
果 。 图 11 本 文 算法 得 到 的 全 球 风 场 可 视 化 效果 
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3.3 ”算法 效率 对 比 实验 


3.3.1 与 其 


取 
简 
为 


也 可 视 化 算法 的 对 比 实验 
为 了 获取 本 文 算法 与 其 他 可 视 化 算法 的 时 耗 差 距 ， 实 验 选 
了 文献 [1, 2] 中 所 提 的 点 图 标 法 、 矢 量 线 法 以 及 拓扑 法 (分别 
写 为 PIM、VLM、TM) 与 本 文 算法 (NGCLIC) 进行 对 比 实验 。 表 1 
上 述 可 视 化 方法 的 用 时 统计 。 
表 1 不 同 可 视 化 方法 的 用 时 比较 /s 
TM NGCLIC 
720x361 102 136 241 2.03 
1440x481 2.54 345 532 4.39 
从 时 间 上 可 以 看 出 ， 本 文 算法 的 时 耗 比 点 图 标 法 和 矢量 线 
比 拓 扑 法 要 低 ， 但 在 本 文 算法 得 到 的 可 视 化 效果 中 ， 纹 


矢量 场 PIM VLM 


上 的 点 之 间 上 共有 更 强 的 相关 性 ， 这 种 特征 不 仅 体现 了 矢 


场 数据 的 连续 性 ， 而 且 达 到 纹理 模拟 真实 场景 的 远 真 效果 ， 
图 标 法 和 矢量 线 法 无 法 达到 的 。 


3.3.2 与 其 他 LIC 改进 算法 的 对 比 实验 


实 
为 
表 


比 
中 
运 
轻 
式 
会 
新 


YE 


加 
因 


4 


射 


法 依旧 存在 时 耗 问 题 。 此 外 , 本 文 算法 没有 结合 GPL 


为 了 分 析 本 文 算法 与 其 他 LIC 改进 算法 在 性 能 上 的 差异 ， 
验 又 选取 了 文献 [5 一 9] 中 所 提 到 的 LIC 改进 算法 (分 别 简写 
HFLIC、HSVLIC、LTLIC、NCLIC) 与 本 文 算法 进行 对 比 实验 。 
2 为 不 同 的 LIC 改进 算法 用 时 统计 。 

表 2 不 同 的 LIC 改进 算法 用 时 比较 /s 

REH HFLIC HSVLIC LTLIC NCLIC NGCLIC 
720x361 2.41 2.23 2.14 1.75 2.03 

1440x481 5.32 4.94 4.68 4.27 4.39 


从 表 2 可 以 看 出 ， 本 文 算法 在 保证 可 视 化 的 质量 的 同时 ， 
文献 [5, 7, 8] 中 提 到 的 算法 在 效率 上 有 一 定 提升 ,但 比 文献 [9] 
的 算法 效率 要 低 。 本 文 算 法 通过 运用 FastLIC 思想 实现 LIC 
算 模块 、 划 分 输入 纹理 为 n 个 模块 同步 执行 LIC 运算 、 选 用 
量 级 的 OpenCV 处 理 引 擎 、 并 行 实现 LIC 运算 和 非 线 性 渐变 
颜色 映射 等 一 系列 措施 来 提高 算法 的 运行 速度 ， 这 种 优势 也 
随 着 数据 量 的 增加 而 慢 慢 显现 ; 但 本 文 的 LIC 模块 需要 生成 
的 纹理 ， 分 模块 同步 计算 需要 系统 开销 ， 且 算法 的 合成 模块 
要 同时 等 到 LIC 结果 和 颜色 映射 结果 才能 执行 ， 所 以 本 文 算 
方式 进行 
速 泻 染 , 这 也 是 本 文 算法 比 文献 [9] 中 提 到 的 算法 效率 低 的 原 
， 但 本 文 算法 对 硬件 要 求 不 高 ， 适 用 性 更 广 。 


结束 语 


本 文 针 对 经 典 LIC 算法 和 LIC 改进 算法 在 矢量 场 强度 上 了 映 
弱 的 缺点 ， 提 出 了 基于 非 线 性 渐变 式 颜色 映射 的 LIC 改进 算 


法 


。 通 过 实验 验证 ， 改 进 后 的 算法 适合 变化 剧烈 的 海洋 流 场 、 


以 


场 等 典型 的 矢量 场 ， 不 仅 同时 展现 了 矢量 场 的 方向 和 强度 ， 


n 
的 


彩 
段 


且 能 够 清晰 、 准 确 地 反映 出 矢量 场 全 方位 信息 以 及 局 部 区 域 
变化 趋势 。 所 绘制 出 的 彩色 纹理 连续 平滑 ， 保 持 了 较 高 的 色 
亮度 和 视觉 分 辨 率 ， 取 得 了 较 好 的 可 视 化 效果 。 用 可 视 化 手 
处 理 及 表现 矢量 环境 已 成 为 了 迫切 需要 解决 的 课题 。 本 文 算 
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法 在 海洋 流 场 、 风 场 可 视 化 中 的 应 用 ， 也 为 气动 流 场 、 磁 极 场 
等 其 他 矢量 场 提供 新 的 研究 方法 。 未 来 本 文 还 将 在 此 基础 上 充 


分 结合 数据 挖掘 的 方法 去 探寻 矢量 场 的 变化 规 得 


和 特征 结构 以 


pu 


及 其 他 现象 ， 并 建立 良好 可 视 化 挖掘 机 制 。 
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